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Abstract — This paper presents an analytical model and a 
geometric numerical integrator for a tethered spacecraft model 
that is composed of two rigid bodies connected by an elastic 
tether. This model includes important dynamic characteristics 
of tethered spacecraft in orbit, namely the nonlinear coupling 
between deformations of tether, rotational dynamics of rigid 
bodies, a reeling mechanism, and orbital dynamics. A geometric 
numerical integrator, referred to as a Lie group variational in- 
tegrator, is developed to numerically preserve the Hamiltonian 
structure of the presented model and its Lie group configuration 
manifold. The structure-preserving properties are particularly 
useful for studying complex dynamics of a tethered spacecraft 
over a long period of time. These properties are illustrated by 
numerical simulations. 

I. Introduction 

Tethered spacecraft are composed of multiple satellites in 
orbit, that are connected by a thin, long cable. Numerous 
innovative space missions have been envisaged, such as 
propulsion by momentum exchange, extracting energy from 
the Earth's magnetic field, satellite de-orbiting, or Mars 
exploration [1], [2], [3], and several actual missions, such 
as TSS, SEDS, or YES2 by NASA and ESA [1], [4]. 

The dynamics of tethered spacecraft involves nonlinear 
coupling effects between several dynamic modes evolving on 
multiple length and time scales. For example, the length of 
tether typically varies from 20 km to 100 km, but the orbital 
radius of tethered spacecraft is several thousand kilometers. 
The natural frequency of the tether is much higher compared 
to the rotational attitude dynamics or the orbital period of the 
spacecraft. The rotational dynamics of spacecraft is nontriv- 
ially coupled to the tension of the tether, which is affected 
by the reeling mechanism and orbital maneuver. Therefore, 
it is important to accurately model tether dynamics, attitude 
dynamics of spacecraft, reeling mechanisms, gravitational 
force and the interaction between them. 

Several analytic and numerical models have been devel- 
oped for tethered spacecraft. However, due to the com- 
plexities of tethered spacecraft, it is common practice to 
use simplified models. Two point masses connected by a 
rigid tether is considered in [5]. A massless, flexible tether 

Taeyoung Lee, Mechanical and Aerospace Engineering, Florida Institute 
of Technology, Melbourne, FL 39201 taeyoung@f it . edu 

Melvin Leok, Mathematics, University of California at San Diego, La 
Jolla, CA 92093 mleok@math.ucsd.edu 

N. Harris McClamroch, Aerospace Engineering, University of Michigan, 
Ann Arbor, MI 48109 nhm@umich.edu 

*This research has been supported in part by NSF under grants CMMI- 
1029551. 

tThis research has been supported in part by NSF under grants DMS- 
0714223, DMS-0726263, DMS-0747659. 



dynamics is included in [6]. Transverse vibrations of two 
point masses connected by a flexible, but inextensible, tether 
are studied in [7]. These simplified models allow for rigorous 
mathematical analysis, but they may fail to predict the behav- 
iors of an actual tethered spacecraft accurately, particularly 
given the fact that tethered spacecraft operations are based 
on weak nonlinear effects over a long time period. Re- 
cent numerical studies consider more sophisticated tethered 
spacecraft models including a varying tether length. But, in 
these advanced models, rigid body dynamics is ignored [8], 
[9], and a reeling mechanism is neglected [8], [10]. 

The goal of this paper is to develop a high-fidelity analyti- 
cal model and numerical simulations for tethered spacecraft. 
This is an extension of preliminary work that studies a 
string pendulum model with a reeling mechanism [11], 
[12]. The first part of this paper provides a realistic and 
accurate analytical tethered spacecraft model including tether 
deformations, attitude dynamics of rigid bodies, and a reeling 
mechanism. We show that the governing equations of motion 
can be developed using Hamilton's principle. 

The second part of this paper deals with a geometric 
numerical integrator for tethered spacecraft. Geometric nu- 
merical integration is concerned with developing numerical 
integrators that preserve geometric features of a system, such 
as invariants, symmetry, and reversibility [13]. A geometric 
numerical integrator, referred to as a Lie group variational 
integrator, has been developed for a Hamiltonian system on 
an arbitrary Lie group in [14]. 

A tethered spacecraft is a Hamiltonian system, and its 
configuration manifold is expressed as the product of the 
Lie groups SO(3), SE(3), and the space of connected curve 
segments on R 3 . This paper develops a Lie group variational 
integrator for tethered spacecraft based on the results pre- 
sented in [14]. The proposed geometric numerical integrator 
preserves symplecticity and momentum maps, and exhibits 
desirable energy conservation properties. It also respects the 
Lie group structure of the configuration manifold, and avoids 
the singularities and computational complexities associated 
with the use of local coordinates. It can be used to study non- 
local, large amplitude and deformation maneuvers of tethered 
spacecraft accurately over a long time period. 

II. Tethered Spacecraft 

We consider two rigid spacecraft connected by an elastic 
tether. We assume that rigid spacecraft can freely translate 
and rotate in a three-dimensional space, and the tether is 
extensible and flexible. The bending stiffness of the tether is 
not considered as the diameter of the tether is assumed to 
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(a) Reference configuration (b) Deformed configuration 

Fig. 1. Tethered spacecraft model 



be negligible compared to its length. The tether is connected 
to a reeling drum in a base spacecraft, and the other end of 
the tether is connected to a sub-spacecraft. The point where 
the tether is attached to the spacecraft is displaced from the 
center of mass so that the dynamics of the spacecraft is 
coupled to the tether deformations and displacements. This 
model is illustrated in Fig. [T] 



We choose a global reference frame and two body-fixed 
frames. The global reference frame is located at the center 
of the Earth. The first body fixed frame is located at the 
center of mass of the base spacecraft, and the second the 
body-fixed frame is located at the end of the tether where 
the tether is attached to the sub-spacecraft. Since the tether 
is extensible, we need to distinguish between the arc length 
for the stretched deformed configuration and the arc length 
for the unstretched reference configuration. Define 
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the mass of the base spacecraft 

the inertia matrix of the base spacecraft 

the rotation matrix from the first body 

fixed frame to the reference frame 

the angular velocity of the base spacecraft 

represented in its body fixed frame 

the location of the center of mass of the 

base spacecraft represented in the global 

reference frame 

the radius of the reeling drum 

the length of the guideway 

the vector from the center of mass of the 

base spacecraft to the beginning of the 

guideway represented in its body fixed 

frame, p = [d,0,b]. 

the mass of the reeling drum 

the inertia matrix of the reeling drum, 

J r = K r d 2 for a matrix n r G R 3x3 

the total unstretched length of the tether 

the unstretched arc length of the tether 

between the point at which the tether is 

attached to the reeling drum an a material 

point P on the tether 

the stretched arc length of the tether to the 

material point located at s 

the arc length of the tether between the 

point at which the tether is attached to 

the reeling drum and the beginning of the 

guide way 

the deformed location of a material point 

P from the origin of the global reference 

frame; r(s p ,t) = x(t) + R(t)p 

9 = (s p — b — s) jd for s £ [0, s p — b] 

The mass of the tether per unit unstretched 

length 

the mass of the sub-spacecraft 
the inertia matrix of the sub-spacecraft 
the rotation matrix from the second body 
fixed frame to the global reference frame 
the angular velocity of the sub-spacecraft 
represented in its body fixed frame 
the vector from the point where the tether 
is attached to the sub-spacecraft to the 
center of mass of the sub-spacecraft rep- 
resented in its body fixed frame 
control moment applied at the reeling 
drum 



A configuration of this system can be described by the 
locations of all the material points of the tether, r(s,t) 
for s E [0,-L], the location of the base spacecraft, the 
attitude of both spacecraft, and the length of the deployed 
portion of the tether. So, the configuration manifold is G = 
C°°([0,l},R 3 ) x SE(3) x SO(3) x E, where C°°([0, 1], M 3 ) 
denotes the space of smooth connected curve segments on 
M 3 , SO(3) = {R G R 3x3 \R T R = J,det[i?] = 1}, and 
SE(3) = M 3 ©SO(3) [15]. 

Throughout this paper, we assume that: (i) the radius of a 



reeling drum and the length of a guideway is small compared 
to the length of a tether; (ii) the reeling drum rotates about 
the second axis of the first body fixed frame attached to 
the base spacecraft; (iii) the deployed portion of the tether 
is extensible, but the portion of the tether on the reel and 
the guideway inside of the base spacecraft is inextensible; 
(iv) the gravity is uniform over the base spacecraft and the 
sub-spacecraft. 

III. Continuous-time Analytical Model 

In this section, we develop continuous -time equations of 
motion for a tethered spacecraft using Hamilton's variational 
principle. The attitude kinematics equation of the base space- 
craft and the sub-spacecraft is given by 



R = RQ, R s 



Rn s , 



(i) 



where the hat map ■ : M 3 — > so (3) is defined by the condition 
that xy — x x y for any x, y £ M 3 . 

A. Lagrangian 

Kinetic energy: The kinetic energy of the base space- 
craft excluding the reeling drum is given by 



Tf, 1 = -mi • x + — f2 • Jil. 



(2) 



Under the assumption that the radius of a reeling drum is 
much less than the length of the tether, the kinetic energy of 
the reeling drum and the part of the tether inside of the base 
spacecraft can be approximated by 
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T b 2 = ^( m r + LlSp)± ■ X + -pS p S p + -K 2 



(3) 



where k 2 — e 2 • n r e 2 . Let r(s,i) be the partial derivative of 
r(s,t) with respect to t. The kinetic energy of the deployed 
portion of the tether is given by 



T, 



-p.r(s) ■ f(s) ds. 



(4) 



Let p £ M 3 be the vector from the end of the tether to a 
mass element of the sub-spacecraft represented with respect 
to its body fixed frame. The location of the mass element in 
the global reference frame is given by r(L) + R s p. Then, 
the kinetic energy of the sub-spacecraft is given by 
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r(L) + R s {l s p\\ 2 dm 



\m s r(L) ■ r(L) + m s r(L) ■ R s Cl s p s + ^-tt s ■ J S Q S 



(5) 



Here, we use the fact that J B pdm = p c , and J s = 
— J„ p\dra. The total kinetic energy is given by T = 

T bl +T b2 +T t +T s . 



Potential energy: By the assumption that the size of 
reeling drum is small compared to the length of the tether, the 
gravitational potential of the base spacecraft and the reeling 
mechanism is approximated as follows 



Vb = — (m + m r + fj,s p ) 



GM 



(6) 



where the gravitational constant and the mass of the Earth 
are denoted by G and M, respectively. 

The strain of the tether at a material point located at r(s) 
is given by 

e = lira ; = s (s) — 1, 

As-s-O As 

where ( )' denote the partial derivative with respect to s. The 
tangent vector at the material point is given by 
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Since this tangent vector has unit length, we have s'(s) = 
\\r'(s)\\. Therefore, the strain can be written as e = ||r'(s)|| — 
1. Using this, the elastic potential and the gravitational 
potential of the deployed portion of the tether is given by 



V t = 
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where E and A denote the Young's modulus and the sectional 
area of the tether. 

The location of the center of mass of the sub-spacecraft is 
r(L) +R s p s in the global reference frame. Since we assume 
that the gravity is uniform over each spacecraft body, the 
gravitational potential energy of the sub-spacecraft is given 
by 



V, = -TO., 



GM 



iL) + R sPs 



(8) 



From |2])-(|8]l, the Lagrangian of the tethered spacecraft is 
given by 



L = Ti 



bi 



■Tt + Tg — Vb — Vt — V s . 



(9) 



B. Euler-Lagrange Equations 

Let the action integral be © = f t * Ldt. According to 
Hamilton's principle, the variation of the action integral is 
equal to the negative of the virtual work for fixed boundary 
conditions, which yields Euler-Lagrange equations. For the 
given tethered spacecraft model, this requires the following 
three careful consideration: (i) the domain of the integral 
depends on the variable s p (t) at Q; (ii) the rotation matrices 
R, R s that represents the attitudes lie in the nonlinear Lie 
group SO (3); (iii) as the tether is assumed to be inextensible 
in the guideway, and it is extensible outsize of the guideway, 
there exists a discontinuity in strain at the beginning of the 
guideway. 



a) Time-Varying Domain: Due to Q, the variation discontinuity, the speed of the tether changes instantaneously 



of the action integral 5® includes the following term, 
Ito Is M^OO ' ^(X) dsdt. Here, we cannot apply integration 
by parts with respect to s, since the order of the integrals 
cannot be interchanged due to the time dependency in the 
variable s p (t). Instead, we use Green's theorem [11], 

r(s) ■ 6r(s) ds = [ ' [ -r(r(s) • 5r{s)) dsdt, (10) 



p(t) 



where § c represents the counterclockwise line integral on the 



boundary C of the region [to,tj 



x is, 



,(t),L]. The boundary 



C is composed of four lines: (i = to, s € [s p (to),L]), 
(t = t f ,s e [s p (t f ),L]), (t e [t ,t f ],s = s p (t)), and 
(t e [to,tf],s = L). For the first two lines, 5r(s) = 
since t = t ,tf. For the last line, ds — since s is fixed. 
Thus, parameterizing the third line by t, we obtain 



r(s) ■ Sr(s) ds = / r(s p (t)) ■ Sr(s p (t)) s p (t) dt. 

>C J t 

Substituting this into ( fT0| ) and rearranging, we obtain 

r(s) ■ Sf(s) dsdt 



In 



-r(s) ■ Sr(s) ds + r(s p ) ■ Sr(s p ) s p 



dt. 



(ID 



b) Variation of Rotation Matrices: The attitudes of 
spacecraft are represented by the rotation matrix R, R s G 
SO (3). Therefore, the variation of the rotation matrix should 
be consistent with the geometry of the special orthogonal 
group. In [14], it is expressed in terms of the exponential 
map as 



SR=4- 
de 



R € = -T- 



Rexpefj — Rfi, (12) 



for r\ € 



The key idea is expressing the variation of 



a Lie group element in terms of a Lie algebra element. 
This is desirable since the Lie algebra so (3) of the special 
orthogonal group, represented by 3 x 3 skew-symmetric 
matrices, is isomorphic as a Lie algebra to K 3 . As a result, 
the variation of the three-dimenstional rotation matrix R is 
expressed in terms of a vector r] G M 3 . We can directly show 
that {l2} satisfies S(R T R) = SR T R + R T SR = -fj + fj = 
0. The corresponding variation of the angular velocity is 
obtained from the kinematics equation <|TJ : 

d 



sn 



de 
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(13) 



c) Variational Principle with Discontinuity: Let r(s~), 
and r(Sp ) be the material point of the tether just inside the 
guide way, and the material point just outside the guide way, 
respectively. Since the tether is inextensible inside the guide 
way, ||r'(s~)|| = 1. Since the tether is extensible outside 
the guide way, ||r'(s+)|| = 1 + e + , where e + represents the 
strain of the tether just outside the guide way. Due to this 



by the amount e + |s p | at the guide way. 

As a result, the variation of the action integral is not equal 
to the negative of the virtual work done by the external 
control moment u at the reeling drum. Instead, an additional 
term Q, referred to as Carnot energy loss term should be 
introduced [9], [16]. The resulting variational principle is 
given by 



5®+ (Q + u/d)Ss p -ue 2 -'ndt = 0. 



(14) 



to 



The corresponding time rate of change of the total energy 
is given by E = (Q + u/d)s p , where the first term Qs p 
represents the energy dissipation rate due to the velocity and 
strain discontinuity. The corresponding expression for Q has 
been developed in [11]: 

Q = ~M(l|r-'(4)|| - l) 2 s 2 p - ^(||r'( S +)|| - If. 

(15) 

d) Euler-Lagrange Equations: Using these results, and 



the variational principle with discontinuity ( 14 1, we obtain 



the Euler-Lagrange equations for the given tethered space- 
In 



19 1, we require that r(s p ) 



craft model in (16l-(21 
x + Rp for the continuity of the tether. 

These can be simplified in a number of special cases. For 
example, we can substitute s p = when the length of the 
deployed portion of the tether is fixed, and we can set p — 
p s = when the main spacecraft and the sub-spacecraft are 
modeled as point masses instead of rigid bodies. 

IV. Lie Group Variational Integrator 

The Euler-Lagrange equations developed in the previous 
section provide an analytical model for a tethered spacecraft. 
However, the standard finite difference approximations or 
finite element approximations of those equations using a 
general purpose numerical integrator may not preserve the 
geometric properties of the system accurately [13]. 

Lie group variational integrators provide a system- 
atic method of developing geometric numerical integra- 
tors for Lagrangian/Hamiltonian systems evolving on a Lie 
group [14]. As they are derived from a discrete analogue 
of Hamilton's principle, they preserve symplecticity and the 
momentum map, and it exhibits good total energy behavior. 
They also preserve the Lie group structure as they update a 
group element using the group operation. These properties 
are critical for accurate and efficient simulations of complex 
dynamics of multibody systems [17]. 

In this section, we develop a Lie group variational inte- 
grator for a tethered spacecraft. 

A. Discretized Tethered Spacecraft Model 

Let h > be a fixed time-step. The value of variables 
at t = t + kh is denoted by a subscript k. We discretize 
the deployed portion of the tether using N identical line 
elements. Since the unstretched length of the deployed 
portion of the tether is L — s Pk , the unstretched length of 



each element is h 



L-i 



N 



-. Let the subscript a denote the 
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where F(s) = EA 1 ,. ,,_!! ,. r'(s) denotes the tension of the tether. 
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variables related to the a-th element. The natural coordinate 
on the a-th element is defined by 



Cfc.aO) 



l)h 



(22) 



for s € [s pt + (a— 1)/^, s Pfc + a/fc]. This varies between and 
1 on the a-th element. Let Sq , Si be shape functions given 
by Sb(C) — 1 — C> an d <Si (C) — C These shape functions are 
also referred to as tent functions. 

Using this finite element model, the position vector r(s, t) 
of a material point in the a-th element is approximated as 
follows: 



r k (s) = S (( k ,a)rk,a + S 1 (Ck,a)rk.a 



+1- 



(23) 



Therefore, a configuration of the presented discretized 
tethered spacecraft at t = kh + to is described by g k — 
(xk; R k ; s Pk ;r k ,i, . . . ,r k ,N+r, R 8k ), and the corresponding 
configuration manifold is G = R 3 x SO(3) x R x (M 3 ) JV+1 x 
SO (3). This is a Lie group where the group acts on itself by 
the diagonal action [15]: the group action on x k , s Pk , and 
r k , a is addition, and the group action on R k ,R Sk is matrix 
multiplication. 

We define a discrete-time kinematics equation 
using the group action as follows. Define f k = 

(Ax k ;F k ;As Pk ;Ar kA ,..., Ar k>N+ f, F Sh ) G G such 
that g k+1 = g k f k : 

(x k+ i; R k+ i; s Pk+1 ;r k+ i, a ; R k+ i) = 

(x k + Ax k ;R k F k ; s Pk + As Pk ;r k . a + Ar k , a ;R Sk F Sk ). 

(24) 

Therefore, f k £ G represents the relative update between two 
integration steps. This ensures that the structure of the Lie 
group configuration manifold is numerically preserved since 
g k is updated by f k using the right Lie group action of G 
on itself. 

B. Discrete Lagrangian 

A discrete Lagrangian Ld(g k ,f k ) ■ G x G — s- K is an 
approximation of the Jacobi solution of the Hamilton-Jacobi 
equation, which is given by the integral of the Lagrangian 



along the exact solution of the Euler-Lagrange equations over 
a single time-step: 

r-h 



Ld{g k ,f k ) 



L(g(t),g-\t)g(t))dt, 



where g(t) : [0, h] — > G satisfies Euler-Lagrange equations 
with boundary conditions <?(0) = g k , g(h) = g k f k . The 
resulting discrete-time Lagrangian system, referred to as 
a variational integrator, approximates the Euler-Lagrange 
equations to the same order of accuracy as the discrete 
Lagrangian approximates the Jacobi solution [18]. 

We construct a discrete Lagrangian for the tethered space- 
craft using the trapezoidal rule. From the attitude kinetics 
equations the angular velocity is approximated by 



1 rp 1 

R k (Rk+i - Rk) = T\F k 



ft KV w h 

Define a non-standard inertia matrix Jj = 



I)- 



|tr[J]J— J. Using 
the trace operation and the non-standard inertia matrix, the 
rotational kinetic energy of the main spacecraft at Q can 
be written in terms of ft as gf2 • JQ = itr[f2Jd" ]■ Using 
this, and ([3]), the kinetic energy of the base spacecraft is 
given by 



T, 
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(m + m r + fis Pk )Ax k ■ Ax k 



2h 2 



h^ 1 



F k )J d ], (25) 



where we use properties of the trace operator: tr[AB] = 
tv[BA] = ty[A T B T ] for any matrices A, B e M 3x3 . 

Next, we find the kinetic energy of the tether. Using the 



chain rule, the partial derivative of r k (s) given by (23 1 with 
respect to t is given by 

r k (s) = ^-j 5 (Cfc,a)Ar fc!tI + 5i(Cfc,a)Ar fc!a+ i 



(L-s) {r k>a - r k , a+ i) 



(L 
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Substituting this into Q, the contribution of the a-th tether 
element to the kinetic energy of the tether is given by 
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MlAr k . a+1 ■ Ar k ,a+l 
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M^ a As Pk 
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(26) 



where inertia matrices are given by 

Ml = \jilk, Ml=Ml 

1 (37V 2 + 3A^ + l-67Va-3a + 3a 2 ) 



M 3i _ l^(2 + :! V-> 



1U k,a 



iV 2 

l_(l + 3AT-3a) 
6™ Af 

(rk,a - rk, a +i)- 



6^ AT 

Similar to (25 1, from (|5j, the kinetic energy of the sub- 
spacecraft is given by 



T k ,. 
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Y 2 



m s Ar k ^ N+1 ■ R Sk (F Sk - I)p s 



(27) 



From ( 25 I, ( 26 1, ( 27 1, the total kinetic energy of the dis- 
cretized tethered spacecraft is given by 



n = T, 



k.b 



N 

£ 

a=l 



Tk,a + T k 



(28) 



Similarly, from ([6]), (|7j, and (jSJ, the total potential energy 
is given by 

-GM(m + m r + ~ps Pk ) - 



V k 



\Xk\ 
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1 iL4 

2 l k 
GMm, 



—2GM~pl 



1 



||ffc, a + rfc, a+ l|| 

rk,a+i - r k , a \\ - h) 2 

1 

lkfe,iv+i + R Sk p s \\ ' 



(29) 



Using ( 28 1, ( 29 1, we choose the discrete-Lagrangian of the 
discretized tethered spacecraft as follows: 

L d k (9k,fk) = hT k (gk,fk) - 7jV k (g k ,fk) - ^V k+ i(g k , f k ). 

(30) 

C. Discrete-time Euler-Lagrange Equations 

We define the discrete action sum 0^ = 
X]fe=i Ld k (<7fe, fk)- According to the discrete Hamilton's 
principle, the variation of the action sum is equal to the 
negative of the discrete virtual work. This yields discrete- 
time Euler-Lagrange equations, referred to as variational 
integrators. 

In [14], the following Lie group variational integrator has 
been developed for Lagrangian systems on an arbitrary Lie 
group: 

Ad*-, • (T:L A • T> fh L dk ) 

U dk + Q dk = 0, 



(3D 



9k+l ~ 9k fk, 



(32) 



where TL : TG — > TG is the tangent map of the left 
translation, Df represents the derivative with respect to /, 
and Ad* : G x g* — > g* is co-Ad operator [15]. 

The virtual work due to the control input and the Camot 
energy loss are denoted by U dk G 0*, and Q dk E g*, 
respectively, and they are chosen as 

-u k Ss k - hu k e 2 ■ rjk, 



U dk ■ (g-'Sgk) 
Qd k ■ (g^Sgk) 



(jlAstJh'+EA) 



|r-fc,2 - r fe) i|| - l k ) 5s 



(33) 



(34) 



By substituting |30|, ([33j, and (|34| into ([3TJ and we 
obtain a Lie group variational integrator for the given discrete 
tethered spacecraft model. Due to page limits, we do not 
present the detailed results in this paper. But, we demonstrate 
their computational properties in the next section. 

V. Numerical Example 
Properties of the tethered spacecraft are chosen as 

m = 490 kg, m r = 10 kg, m s = 150 kg, 
I = 120 km, p = 24.7 kg/km, EA = 659700 N, 
J = diag[5675.8, 5675.8, 6125] kgm 2 , p = [0.5, 0.0, l]m, 
J s = diag[500, 500, 300] kgm 2 , p s = [0, 0, 1] m. 

Initially, the base spacecraft is on a circular orbit with an 
altitude of 300 km, and the tether and the sub-spacecraft 
are aligned along the radial direction. The initial unstretched 
length of the deployed portion of the tether is 20 km, i.e., 
s p (0) — 100 km. The initial velocity at each point of 
the tether and the sub-spacecraft is chosen such that it 
corresponds to the velocity of a circular orbit at their altitude. 

We consider the following three cases. In the first case, 
the reeling drum is fixed so that the length of the deployed 
portion of the tether is fixed, i.e., s p = 0. In the second case, 
the reeling drum is free to rotate, and the tether is released by 
gravity to 100 km. The third case is the same as the first case 
except that the initial velocities of the base spacecraft and 
the sub-spacecraft are perturbed by about 15% to generate a 
tumbling motion. These are summarized with the time-step 
and the simulation time as follows: 
Description 



h(s) 


t f (s) 


0.05 


6000 


0.05 


3848 


0.01 


500 



Case 1 Fixed reeling drum 
Case 2 Releasing the tether to 100 km 
Case 3 Velocity perturbation of Case 1 

Note that the orbital period of a point mass on the circular 
orbit with the altitude of 300 km is 5410 seconds. For all 
cases, the number of tether elements is N = 20. 

The following figures illustrate simulation results for each 
case. We consider a fictitious local vertical, local horizontal 
(LVLH) frame that is attached to an imaginary spacecraft on 
a circular orbit with an altitude of 300 km. For each figure, 
we have the following subfigures: (a) the maneuvers of the 
tethered spacecraft are illustrated with respect to the LVLH 
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VI. Conclusions 

We develop continuous-time equations of motion and a ge- 
ometric numerical integrator for a tethered spacecraft model 
that includes tether deformation, spacecraft attitude dynam- 
ics, and a reeling mechanism. This provides an analytical 
model that is defined globally on the Lie group configuration 
manifold, and the Lie group variational integrator preserves 
the underlying geometric features, thereby yielding a reliable 
numerical simulation tool for complex maneuvers over a long 
time period. 
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frame. To represent the attitude dynamics of spacecraft, the 
size of the spacecraft is increased by a factor of 100, and 
the relative strain distribution of the tether at each instant 
is represented by a color shading (animations illustrating 
these maneuvers are also available at |http: / /my. fit 7] 
|edu/ -taeyoungj l. The remaining subfigures show: (b) the 
energy transfer, (c) the computed total energy deviation 
from its initial value, (d) the angular velocity of the base 
spacecraft, and (e) the unstretched/stretched length of the 
tether. 

In the first case, we observe a pendulum-like motion 
where the tether is taut and its stretch length is almost 
close to the unstretched length. But, there exists a strain 
wave that propagates along the tether, and nontrivial attitude 
maneuvers for the base spacecraft and the sub-spacecraft. 
The proposed Lie group variational integrator exhibit ex- 
cellent conservation properties: the maximum relative total 
energy deviation is 2.37 x 10~ 8 % of its initial value, and 
the maximum orthogonality error of rotation matrices is 
max{||/ - R T R\\} = 1.32 x 10~ 13 . 

In the second case, the tether is deployed by gravity 
gradient effects, and due to the Carnot energy term discussed 
in the previous section, the total energy increases slightly. As 
the mass in the base spacecraft is transferred to the deployed 
portion of the tether, there is a transfer of kinetic energy 



between two parts, as seen in Fig. 3(b) 



The third case is most challenging: there are in-plane 
and out-of-plane tumbling maneuvers, while the tether is 
stretched by 25 %, and the attitude dynamics of spacecraft 
is nontrivially excited with a large angular velocity. The 
proposed Lie group variational integrator computes the com- 
plex dynamics of this tethered spacecraft accurately. The 
maximum relative total energy deviation is 3.48 x 10" 4 %, 
and the maximum orthogonality error of rotation matrices is 
max{|U - R T R\\} = 8.03 x 10~ 14 . 





(a) Snapshots observed at the LVLH frame (km) (The size of spacecraft 
is increased by a factor of 100 to illustrate attitude dynamics.) 



(a) Snapshots observed at the LVLH frame (km) (The size of spacecraft 
is increased by a factor of 100 to illustrate attitude dynamics.) 
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Fig. 2. Case 1: Circular orbit, Fixed unstretched tether length 



Fig. 3. Case 2: Circular orbit, Releasing tether 
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(a) Snapshots observed at the LVLH frame (km) (The size of space- 
craft is increased by a factor of 100 to illustrate attitude dynamics.) 
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(d) Angular velocity of the base space- (e) Unstretched length of the deployed 
craft Q part of the tether (red), stretched 

length (blue) 



Fig. 4. Case 3: Perturbed circular orbit, Fixed unstretched tether length 



